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ABSTRACT 

We examine the r-process in the neutrino-driven proto-neutron-star (PNS) wind of core-collapse 
supernovae in light of the recent findings of massive neutron stars in binaries as well as of an indication 
of neutron-richness in the PNS ejecta because of the nucleon potential corrections on neutrino opacities. 
To this end, a spherically symmetric, general relativistic, steady-state wind model is applied for a 
wide range of PNS masses between 1.2M Q and 2.4M Q with the latter reaching the causality limit. 
Nucleosynthesis calculations with these PNS models are performed by assuming a time evolution of 
electron fraction with its minimal value of Y c = 0.4, which mimics recent hydrodynamical results. 
The fundamental nucleosynthetic aspect of the PNS wind is found to be the production of Sr, Y and 
Zr in quasi-equilibrium and of the elements with A w 90-110 by a weak r-process, which can be an 
explanation for the abundance signatures in r-process-poor Galactic halo stars. PNSs more massive 
than 2.0Mq can eject heavy r-process elements, however, with substantially smaller amount than 
what is needed to account for the solar content. PNS winds can be thus the major origin of light 
trans-iron elements but no more than 10% of those heavier than A ~ 110, although they may be the 
sources of the low-level abundances of Sr and Ba found in numerous metal-poor stars if the maximum 
mass of PNSs exceeds 2.0M©. 
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1. INTRODUCTION 

Proto-neutron-star (PNS) wind of core-collapse super- 
novae (CCSNe), the outflows driven by neutrino heat- 
ing, has long been suggested to be the major site of 
the r-p rocess (rapid neutron-capture process) s ince early 
1990's (IMever et al.l[l99l IWoosiev et al.lU99l . One of 
the problems in these early works was the very high 
entropy in the wind, S ~ 400 fee nucleon -1 (&b is the 
Boltzmann constant), which was not confirmed by sub- 
sequent works (S ~ 100 fen n ucleon" , ITakahashi et al.l 
119941: iQian fc Wooslevl 11996( 1. Gener al relativistic ef- 
fects were found to increase entropy (jCardall fc Fuller! 
fl99l but needed a PNS more massive tha n M ~ 2.0Mp, 
for robust r-processing (lOtsuki et al.l2000HWanaio et alj 
120011: IThompson et al.l 12001( 1. Another problem was 
an unacceptable ov erproduction of som e species such 
as Sr , Y, and Zr (jWooslev et al.1 119941 : IWanaio et all 
2001). More seriously, hydrodynamical simulations 
of CCSNe with elaborate neutrino transport indicated 
proton- richness in th e wind ejecta ((Fischer et al.l l2010t 
iHiidepohl et al.ll2010H . These works seemed to exclude 
the PNS wind scenario as the r-process site. 

Recent works on the effect of nucleon potential cor- 
rections for neutrino opaci ties seem, in part, to re- 
vive t he PNS wind scenario (iReddv et al.l 1998: R oberts! 



2012t iMartmez-Pinedo et al.l l2012t IRoberts et all 120121: 



Horowitz et al.ll2012f) . These works predict that the elec 



tron fraction (Y e ; number of protons per nucleon) drops 
off from an initially proton-rich value to the minimal 
value of ~ 0.42-0.45 and increases again towards ~ 0.5 
in the late wind phase, the behavior not expected in pre- 
vious works. Recent discoveries of massive NSs in binary 
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systems with a precision me asurement of M = 1. 97 ± 
0.04M Q for PSR J1614-2230 (IDemorest et aLllMtCh and 
an inferred mass of M ~ 2.4M Q for PSR B1957+20 
(van Kerkwiik et aDl20ll are also encouraging for the 
PNS wind scenario. 

In this Letter we aim to revisit the issue of the r-process 
in PNS winds in light of these recent findings. This is to 
extend and improve the previous nucleosynthesis studies, 
which were based on limited hydrodynamical outco mes 
(e.g., IWoosiev etHl Il994t lArcones fc Montesl f201lT ) or 
semi-analytic solutions with a few selected parameter 
sets (e.g., M = 1.4M and 2.0M ? : IWanaio et al.ll200l 
with Y c evolutions that were incompatible with the re- 
cent works. To this end, a semi-analytical wind model 
(|Wanaio et al.l I2001T 1 is applied for a wide range of M 
between 1.2M Q and 2.4M (Section 2). Nucleosynthe- 
sis calculations are performed with the wind solutions by 
assuming a ti me evolution o f Y e (S ection 3) that mimics 
the result of IRoberts et al.l (|2012f l ■ The nucleosynthe- 
sis yields are mass-integrated to compare with the solar 
r-process abundances as well as those in Galactic halo 
stars. We then discuss whether PNS winds can be the 
sources of r-process elements in the Galaxy. 

2. WIND MODEL 

Although the underlying physics what causes the 
explosions of CC SNe has been under debate (e.g., 
Uanka et al.l I2012T) . it is known that the neutrino-driven 
outflows after evacuation of the early convective ejecta 
are well described by the steady-state (semi-) ana- 
lytica l solutions of PNS wind models (|Duncan et al.l 
1981 IQian fc Wooslevl ll9M ICardall fc Fuller! 11997 



Otsuki et al.H2000t IWanaio . et al.ll200U IThompson et alj 
20011) . In this study, we use the spherically symmet- 
ric, general relativ istic, semi-analytic wind model in 
IWanaio et al.l (|200lL for more detail, see their Section 2). 
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Fig. 1. — Time evolutions of L„ (top left) and R (top right) adopted in this study. Resulting S (middle left), r (middle right), S/r 1 / 3 
(bottom left), and M (bottom right) are shown as functions of t. In the bottom- left panel, the Y e 's, above which the production of A ~ 200 
nuclei are predicted, are indicated by dashed lines. 



The average neutrino energies are taken to be 12, 14, and 
14 MeV for electron neutrino, electron antincutrino, and 
heavy lepton neutrinos, respectively (acc ording to, e.g., 
IFischer et a!1l2010t iHudepohl et alil2010D . The equation 
of state for ions (ideal gas) and arbitrarily degenerate, ar- 
bitrarily relativisti c electr ons and positrons is taken from 
iTimmes fc Swestvl (pOOOt ). 

Each wind (i.e., transonic) solution can be obtained 
for a given set of (M,R,L V ); we assume the PNS ra- 
dius R to be the same as the neutrinosphere and the 
neutrino luminosities of all the flavors to have the same 
value L v . We consider the models of M/Mq = 1.2, 
1.4, 1.6, 1.8, 2.0, 2.2, and 2.4, which co yer the range of 
measured and estimated masses of NSs (|Demorest et al.l 



120101 : Ivan Kerkwijk et aTll20lH lLa"ttimerll2011| ). For L v 

and R, phenomenological time evolutions during the first 
10 s after core bounce (to = 0.20 s < t < t\ = 10 s) 
are adopted as follows. To roughly mimic recent re- 
sults of lon g-term simulations over the PNS coolin g 
phase (e.g.. IFischer et al.l [Mot IHudepohl et al.l 120101 ), 
we assume L v (t) = L„,o(£Ao) _1 with L v $ (erg s _1 ) = 
10 52 - 4 = 3.98 x 10 52 (Fig. H top left). We also assume 
R(L V ) = (Ro - Ri)(L v /L vfi ) + Ri with R = 30 km 
and R\ = 10 km so that each wind solution can be 
obtained from a given set of (M,L V ). This is equiva- 
lent to set R(t) = (R - R^t/to)- 1 + R ± (Fig. Q] top 
right). Note that R(t±) — 10.4 km matches the lower 
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Fig. 2. — Top: time evolution of Y c adopted for nucleosyn- 
thesis calculations (y Cjm i n = 0.40; solid curve). The case with 
Ye min = 0.45 is also shown by the dashed curve. Bottom: values 
of /200 (Eq. [1]) as functions of t. The cases for y Cjm i n = 0.40 
and 0.45 are shown by solid and dashed curves, respectively. The 
solid and dashed horizontal lines indicate /200 = 1 and /130 = 1, 
respectively, above which the production of A ~ 200 and A ~ 130 
nuclei are expected. 

bound of the constraint for cold NSs (with M = 1 . 4M^) , 
10.4km < R < 12.9km, inferred by lSteiner et al.1 (|2013[ i. 
With this assumption, the 2.4M Q model reaches the 
smallest radius (w 10.3 km) allowed by the causality 
limit (the speed of sound mu st not exceed t he speed of 
light), R > 4.3 (M/Mq) km (|Lattimerll20llT) . The PNS 
with M = 2.4M Q should be thus taken as the absolute 
extreme model0 

Wind solutions for each M model are computed for 
log (erg s" 1 ) = 52.60, 52.59, 50.90 (t < t < h; 
171 L^'s). The middle and bottom panels (Fig. [T]) 
illustrate the resultin g basic properties . We confirm 
the previous results (lOtsuki et al.l 120001 : iWanaio et al.l 
[200lt iThompson et al.l 1200 ID that the asymptotic en- 
tropy (S; middle left) increases with time, being sys- 
tematically greater for more massive PNSs. We further 
find a strong sensitivity of S to M for > 2.0M Q , which 
reaches 338 k B nucleon" 1 for M = 2.4M . This is a con- 
sequence of the general relativistic effects that are par- 
ticula rly important when M IR is close to the causality 
limit flCardall & Fullerl [19971) . We also find systemati- 

2 IWa naio ct al. (2001) showed that, with inclusion of general 
relativistic effects, nucleosynthetic outcomes were roughly scaled 
with M/R. Nucleosynthetic results of M/Mq = 1.2-2.4 with Ri = 
10 km here would be thus similar to those of, e.g., M/Mq = 1.4—2.9 
with Ri = 12 km. 



cally smaller expansion timescales (r; defi ned as the e- 
foldin g time of temperature below 0.5 MeV. lOtsuki et all 
21100! for more massive PNSs, which take minimal val- 
ues at t ~ 1-2 s and increase with time (Fig. [T] middle 
right). 

As analytically shown by Hoffma n et al.1 ((1997), the 
ratio S/t 1 / 3 (with a fixed Y c ) serves as the measure of 
the strength of r-processing. We find in the bottom-left 
panel (Fig. Q} that S/t 1 / 3, increases with time and sat- 
urates at t ~ 5 s as a result of the increasing r that 
counterbalances the increasing S. That is, the strength 
of r-processing becomes mostly independent of time for 
t > 5 s if Y c is kept constant. The dashed lines in- 
dicate the values of Y e above which the production of 
A ~ 200 nucl ei are expected , accor ding to the analyti- 
cal formula in iHoffman et al.l (|1997l their Eq. (19)). We 
find that an unacceptably low Y c (< 0.25) is required for 
M = 1.4Mq. If a currently predicted minimal value of 
Y c ~ 0.42-0.45 were taken, one would need a PNS with 
the mass close to the extreme case of M = 2.4M Q for 
the production of heavy r-process nuclei0 The bottom- 
right panel shows the mass ejection rates (M) that are 
systematically smaller for more massive PNS models and 
quickly decrease with time. This indicates that the wind 
ejecta are dominated by the early components with small 
S/t 1 / 3 . The very late ejecta for t > 10 s (if any; not con- 
sidered in this study) would be unimportant. 

The time evolution of Y Cl which is needed for 
nucleosynthesis calculations, is assumed as Y c (t) = 
ciCOsh[c 2 (i - imin)] + c 3 (Fig. [2] top), where c 2 = 1.0 



for t < t. 



3.0 s and c 2 = 0.10 for t > t nnn . The coef- 



ficients ci and C3 are determined to satisfy V c (^min) = 
Y^min and, for t < t m - m and t > t m - ln , respectively, 
Y e (to) = 0.55 and Y e (ti) = 0.50. This roughly m imics the 
hydrodynamical result by iRoberts et al.l ()2012l . see their 
Fig. 5). We adopt Y^ m in = 0.40 (solid curve in the top 
panel of Fig. [2j, the value s lightly smaller than ~ 0.42- 
0.45 in IRoberts efall (pOll ). While Y c drops m response 
to the PNS contraction, the increasing neutrinospheric 
density suppresses the charged current neutrino inter- 
actions by Pauli blocking and Y e cannot decrease at 
late times dFischer et al.ll2012D. Note that the a-effect 
(|McLaughlin et a,l. | ll996t iMever et al.lll998D . which were 
not considered in IRoberts et al.1 (|2012l ). would slightly 
shift Y c towards ~ 0.5. The value of l^ m in here may 
thus be taken as the absolute lower limit for PNS winds. 

The lower panel (Fig. [5]) shows the condition for 
making the th i rd pe ak nuclei (A ~ 200) according to 
IHoffman etaLI (|1997t ). 



(5/230 fc B nucleon _1 ) 
(y c /0.40)(r/20ms) 1 /3 



> 1, 0.38 <Y C < 0.46(1) 



This reflects the value of Y e in addition to the com- 
bination S/t 1 / 3 (Fig. [U bottom left). We find that 
only the extreme model of M = 2.4M Q satisfies this 
condition (the region above the horizontal solid line). 
Also indicated by the horizontal dashed line is the con- 



3 lOtsuki et at] ||2000T) and lWanaio et al.l I T20011 ) showed that the 
wind of (M,L V ) = (2.0M Q , 10 52 erg s _1 ) with Y c = 0.40 led to a 
robust r-process. Their PNS radius was, however, fixed to 10 km, 
which was appreciably smaller than our more reasonable, time- 
dependent value of R(L L , = 10 52 erg s _1 ) = 15 km. 
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TABLE 1 

EjECTA MASSES (iN UNITS OF 1O _5 M0) 



M/Mq 


1.2 


1.4 


1.6 


1.8 


2.0 


2.2 


2.4 


total 


219 


143 


100 


74.1 


56.7 


44.6 


36.0 


4 He 


122 


92.7 


71.9 


56.9 


45.8 


37.4 


31.0 


A > 100 


2.19 


2.75 


2.76 


2.27 


1.78 


1.37 


0.893 


Sr 


3.61 


1.92 


1.09 


0.627 


0.346 


0.177 


0.0764 


13a 


0.00 


0.00 


0.00 


0.00 


0.0420 


0.0373 


0.0199 


Eu 


0.00 


0.00 


0.00 


0.00 


0.00452 


0.00585 


0.00305 



(GT2 [Tachibana et al.l I1990D obtained with the HFB- 
21 masses. Neutrino interactions, which would slightly 
shift Y c by the a-effect, are not included. Using ther- 
modynamic trajectories of PNS winds, the calculations 
are started when the temperature decreases to 10 GK, 
assuming initially free protons and neutrons with mass 
fractions Y c and 1 — Y e , respectivelyO 

The nucleosynthetic abundances are mass-integrated 
(Fig. [3J top) by adopting M for each PNS model. For 
comparison purposes, the solar r-process compositions 
(circles) are also plotted to match the third peak height 
(A ~ 195) for the M = 2.4M© model. As anticipated 
from the lower panel of Figure [5J only the extreme model 
of M = 2AMq satisfactorily accounts for the production 
of heavy r-process nuclei up to Th (A = 232) and U {A = 
235 and 238). The 2.2M model reaches the third peak 
abundances but those beyond. The 2.OM model reaches 
the second (A ~ 130) but the third peak abundances. 
We find no strong r-processing for the models with M < 
2.0M Q . 

We find, however, quite robust abundance patterns 
below A ~ 110, which appears a fundamental aspect 
of nucleosynthesis in PNS winds. The double peaks at 
A « 56 and 90 with a trough between them are formed 
in quasi- nuclear equilibrium (QSE; > 4 GK). Note also 
that the overproduction of iV = 50 species 8 8 Sr, 89 Y, 
90 Zr (IWooslev et alJll99i IWanaio et al.l l200lh are not 
prominent in our result. This is due to the short du- 
ration of moderate S (< 100 fee nucleon -1 ; Fig. [1]) with 
Y c ~ 0.45 (Fig. [31), in which the N = 50 species copiously 
form in QSE. The lower panel of Figure [3] shows the ra- 
tios of nucleosynthetic abundances relative to their solar 
r-process values (normalized at A — 90). For 2.2M and 
2.4M models, the ratios are more or less flat between 
A = 90 and 200, although deviations from unity are seen 
everywhere. 

Table 1 provides the masses (in units of 1O _5 M ) of the 
total ejecta, 4 He, those with A > 100, Sr, Ba, and Eu, for 
all the PNS models. The total ejecta masses span a factor 
of 6 with smaller values for more massive PNSs. The 
larger fractions of 4 He in more massive models, however, 
lead to the ejecta masses for A > 100 (total masses of re- 
process nuclei) ranging only a factor of 2.5. The masses 
of Sr range a factor of 50 with the greater amount for 
less massive models. Ba and Eu are produced only in 
the massive models with M > 2.OM . 

Studies of Galactic chemical evolution estimate the av- 
erage mass of Eu per CCSN event (if they were the ori- 
gin) to be ~ 10" 7 M Q (jlshimaru fc Wanaiol[l999l ). that 
is, ~ a fewlO _5 M for the nuclei with A > 100. Taken at 
the face value, the Eu masses for M > 2.OM reach 30%- 
60% of this requirement. The fraction of events with such 
massive PNSs would be limited to no more than ~ 20% 
of all CCSN events (e.g., > 25M ). The masses of Eu 
from these massive PNSs are, therefore, about 10 times 
smaller than the requirement from Galactic chemical evo- 

4 We examined only the Y e>m j n = 0.40 case with i m i n = 3.0 s. 
Tests showed that a small shift of t m i n did not qualitatively change 
our result. The cases with Y e , m j n = 0.45 corresponded to roughly 
re-scaling M's with ~ 0.1-0.2AfQ smaller values. Note also that 
the presence of th e preceding SN ejec ta that give rise to the 
termination-shocks (Arconcs et al. 2007) do not chang e the gross 
abundance features (Wanaio 2007; Kuroda et al. 2008). 
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Fig. 3. — Top: mass-integrated nuclear abundances, which are 
compared with the solar r-process abundances (circles) that shifted 
to match the third peak height (A ~ 200) for the 2.4M model. 
Bottom: ratios of mass-integrated abundances relative to the solar 
r-process abundances (scaled at A = 90). 

dition for making the second peak (A ~ 130) nuclei, 
/i 30 ~ 1.34/2oo > 1- It indicates that only the mod- 
els with M > 2.0M Q can reach the second peak of the 
r-process abundances. The /200 curves with Yc. m i n re- 
placed by 0.45 are also shown in Figure [2l implying 
slightly weaker r-processing. 

3. NUCLEOSYNTHESIS 

The nucleosynthetic yields for all the (M, L v ) sets 
are computed wit h the reaction network code described 
in IWanaio et al.1 (|2001L I2011bl) . Reaction rates are 
employed from the latest library of REACLIB V2.0 
(|Cvburt et al.l 120101 ) for the experimental evaluations 
when available and th e rest from the theoretical es- 
timates in BRUSLIB dXu et al l I2013D base d on the 
HFB-21 mass predictions (jGorielv et all 120101) . The j3- 
decay rates are taken from the gross theory predictions 
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Fig. 4. — Mass-integrated elementa l abundances compared to the stellar ab undances (top panels) and their ratios (bottom panels). Two 
Galactic halo stars, HD 122563 (left. lHonda et al.|[2006l ; IRoederer et al.ll20"T2T ) and CS 31082-001 (right Igiqueira Mello Jr. et al.|[20l3) are 
taken as representative of r-process-poor/rich stars, respectively. The nucleosynthetic abundances are normalized at Z = 40. 

duces the elements of Z = 46-48 (Pd, Ag, Cd). It could 
be thus possible to interpret that the abundance signa- 
tures of r-process-poor stars were due to a weak r-process 
that reaches Z ~ 50 (M < 2.0M Q ) or 70 [M = 2.OM ) 
with or without additional sources for Z > 50, respec- 
tively. In the right panels, we find that the stellar abun- 
dances between Z = 38 (Sr) and Z = 47 (Ag) are well 
reproduced by massive models with M > 1.6Af . The 
models with M = 2.2Af Q and 2.4M Q produce the heavier 
elements with a similar pattern to that of CS 31082-001 
but with a smaller ratio. Because of the insufficient pro- 
duction of Eu (Table 1), our PNS models would not ac- 
count for the high [Eu/Fe] value in this star. The winds 
from such massive PNSs (M > 2.OM ) could be, how- 
ever, still the source of the low-level abundances (fac- 
tor of several 10 smaller than the avera ge values) of Sr 
and Ba in nume rous metal-poor stars fRoederer 2013; 
lAoki et al.ll2013ft . 



lution (the same holds for Ba). Note that, for massive 
PNS cases, the ejecta masses would b e further reduced 
by fallback or b lack-hole formation (Qian etaUU99a 
iBovd et aLll2012f) . For Sr, the required mass per CCSN 
event is estimated to be ~ 2 x 10~ 6 MfT) from the s olar r- 
process ratio of Sr/Eu = 16.4 (|Sneden et alj |2008). The 
low mass PNS models, which may represent the majority 
of CCSNe, thus overproduce Sr by about a factor of 10. 
The amount of QSE products such as Sr, Y, and Zr is, 
however, highly dependent on the multi-dimensional Y e 
distribution in early times (t < 1 s. lWanaio et al.ll2011aD . 

Figure 2] compares the mass-integrated abundances 
with those of Galactic halo stars. Two well known 
objects are taken as represent ative of r-p r ocess- 
poor (HP 122563, left panels; iHonda et all [20061: 
IRoederer etap 120121) and r-process-rich (CS 31082-001, 
right panels; ISiaueira Mello Jr. et al.l 120131 ) stars with 
the metallicities [Fe/H] = —2.7 and —2.9, respectively. 
These stars have [Eu/Fe] = —0.52 and +1.69, respec- 
tively, well below and above the average value of « +0.5 
at [Fe/H] sa —3. The top and bottom panels show, re- 
spectively, the mass-integrated abundances and their ra- 
tios relative to the stellar abundances, which are normal- 
ized to the stellar abundances at Z — 40. 

In the left panels, we find that the 1.2M Q and 1.4M Q 
models result in reasonable agreement with the stellar 
abundances between Z = 38 (Sr) and Z = 48 (Cd). The 
2.0Mq model nicely reproduces the abundance pattern 
of HD 122563 up to Z = 68 (Er) but somewhat overpro- 



4. CONCLUSION 



We revisited the issue of the r-process in neutrino- 
driven PNS winds in light of recent findings of the mas- 
sive NSs as well as of the neutron-richness in the PNS 
ejecta. Nucleosynthesis calculations were performed with 
the semi-analytical wind models over a wide range of 
the PNS masses (1.2 < M/M < 2.4), assuming a phe- 
nomenological time evolutions of Y c . 

Based on our result, including the extreme model 
(M = 2AMq) that encounters the causality limit, it 
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would be safe to conclude that neutrino-driven PNS 
winds were excluded as the major origin of heavy re- 
process elements. Note that some previous works have 
suggested several mechanisms that help to increase S 
and r educe r, e.g., strong ma gnetic field (jThompsonl 
120031 : iSuzuki fc Nagataki 2 0051) and highly anisotropic 
neutrino emission ( Wanaiol l2006). Our conclusion would 
not be changed if such mechanisms (associated to proba- 
bly rare events) were considered, as far as the dominant 
driving source of the wind is neutrino heating (which 
sets M). Other driving mechanisms such as magneto ro- 
tationally (|Metzger et all 120071; IWinteler et al.ll2012D or 
acoustic-wave driven outflows are beyond the scope of 
this Letter and cannot be excluded as the r-process ori- 



gins on the basis of our result. 

Neutrino-driven PNS winds are, however, promising 
sources that eject light trans-iron elements made in QSE 
(Sr, Y, and Zr) and by a weak r-process (up to Pd, Ag, 
and Cd) , in addition to the early convective ejecta of 
CCSNe (|Wanaio et al.ll2011al ). If not the main origin of 
the r-process elements beyond A ~ 110, they could be 
the sources of low- level abundances of Sr and Ba (e.g., 
[(Sr, Ba)/Fe] < —1) found in extremely metal-poor stars. 
Future surveys or indications of massive NSs, long-term 
hydrodynamical simulations of PNS winds, and galac- 
tic chemical evolution studies will be important to make 
clear the role of PNS winds in the enrichment histories 
of galaxies. 

This work was supported by the JSPS Grants-in-Aid 
for Scientific Research (23224004). 
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